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1. Introduction 

Symbiotic stars are binary stars that are identified spectroscopically as typically having 
three components: a giant, a hot star (usually a white dwarf), and a circumstellar nebula, 
ionized by the hot component's radiation field. Accretion of material from the giant's wind 
onto the white dwarf can lead to eruptions that may last several years and raise the system's 
brightness by several magnitudes. AG Draconis (BD+67°922) is one of the brightest symbi- 
otics at UV and X-ray wavelengths and undergoes outbursts in 3-6 year long intervals that 
repeat every 15 years or so. The observations described here took place during quiescence. 
The components of AG Dra are a KO-4 giant and a hot star, most likely a white dwarf 
(Mikolajewska et al. 1995). For convenience we will refer to the hot component as a white 
dwarf for the remainder of this article. The temperature of the white dwarf is derived through 
fitting a blackbody curve to the continuum observed at X-ray and UV wavelengths, although 
contrasting results are found: Greiner et al. (1997) finding a temperature of 170,000 K from 
ROSAT X-ray spectra, while (Mikolajewska et al. 1995) finding 100,000 K from lUE spectra. 

The system is a non-eclipsing binary with an orbital period of ~550 days (Meinunger 
1979; Galis et al. 1999; Fekel et al. 2000). The radial velocity was found by De Medeiros & 
Mayor (1999) to be —147.46 ± 3.54 km s"-*^ and they give the rotational velocity of the giant 
as 5.9 ± 1.0 km s~^. The abundances of the giant are well determined (Smith et al. 1996) 
and demonstrate that it is metal deficient and possibly a barium star. The distance to AG 
Dra is crucial in fixing many of the properties of the system, yet it is rather uncertain, with 
the Hipparcos satellite only yielding a lower limit of 1 kpc. Mikolajewska et al. (1995) argue 
from a classification of the giant for a distance of 2.5 kpc, which is the value we adopt in the 
present work. 

AG Dra is a particularly attractive target in the far ultraviolet on account of the low 
extinction (Eb-v'=0.05 Mikolajewska et al. 1995) along its sightline, resulting from the 
high latitude of the system {b = 41°). Spectra from the International Ultraviolet Explorer 
{lUE) showed strong emission lines from highly ionized species including Hell A1640, CiV 
AA1548,1550 and Nv A1238,1242 (Viotti et al. 1983). The first reported observations for 
wavelengths below 1200 A were made by Schmid et al. (1999) who identified strong emission 
lines from Ovi A1032,1038 and Hell A1085. 

We report here observations made with the Far Ultraviolet Spectroscopic Explorer (FUSE) 
in the wavelength region 905-1 187 A made in 2000 March and 2001 April during a quiescent 
period of the system. These spectra are at a higher resolution and sensitivity than any 
previous spectra in this range and reveal a number of new plasma diagnostics available for 
interpreting the nebular emission from the system. 
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2. 



The instrument 



FUSE is described in detail in Moos et al. (2000) and Blair & Andersson (2001), and 
we describe the main features of the spectrograph here. There are four telescope channels 
labeled LiFl, LiF2, SiCl, and SiC2, each containing a mirror, a focal plane assembly and a 
grating. The labels refer to the coatings on the mirrors and gratings, either silicon carbide 
(SiC) or lithium fluoride (LiF). The LiF coatings are highly efficient at wavelengths above 
1000 A, but have very little efficiency for smaller wavelengths, hence the SiC channels 
extend the FUSE wavelength range down to 900 A. 

There are two detectors, with the LiFl and SiCl spectra being imaged on detector 1, 
and LiF2 and SiC2 imaged on detector 2. Each detector consists of two segments, labeled 
A and B, which are read out independently by the detector electronics. The LiF and SiC 

spectra on each segment are treated individually by the FUSE calibration pipeline, and so 
one refers to, e.g., the LiFlA spectrum, the SiC2B spectrum, etc. As there is a physical gap 
between each detector's two segments, there is a gap of ~ 10 A between the A and B spectra 
for each channel. 

There are three apertures on the focal plane assembly that are available for observ- 
ing, the largest (30" x 30") and smallest (1.5" x 20") apertures were used for the AG Dra 
observations. 

Guiding is done via a fine error sensor, which is fixed to the LiFl channel. The target 
can move in the remaining channel apertures due to thermally-induced mirror motions. The 
motion of the SiC channels is generally greater than that for LiF2 and occasionally the target 
can leave the aperture. 



Version 1.8.7 of the FUSE calibration pipeline was run on the raw data files. In partic- 
ular, versions 6 and 10 of the fiux and calibration files were used, respectively. For the SiC 
spectra, only that portion of the observation for which the target was in the aperture (see 
Sect. 4) was run through the pipeline to minimize background counts in the spectra. 



The components of AG Dra are not resolved by FUSE and so the spectra presented here 
show emission from the entire system. Two observations were performed through programs 
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Data reduction 



4. 



Observations 
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S312 and P248, and details are given in Table 1. The SiC channels were only ahgned for 
small portions of the total exposure time in 2000 March, leading to effective exposure times 
of 440 s and 300 s for SiCl and SiC2, respectively. The O vi A 1032 flux from the 2000 March 
observation was 2.4 x 10^^*^ erg cm~^ s~^ in the center of the line, exceeding the FUSE 
flux limit of 1.0 X 10^^*^ erg cm~^ s~^ and so a special observing sequence was employed 
for the second observation to ensure the safety of the FUSE detectors. AG Dra was flrst 
observed through the 1.5" x 20" sht (HIRS) which has a throughput around 50% of that 
of the standard 30" x 30" sht and so the count rate on the detector at A1032 was reduced. 
Coupled with the expected fall in the A1032 flux due to observing at a different orbital 
phase of the system (Sect. 13), the effective A1032 flux then fell below the FUSE flux hmit. 
Alignment of the target in the aperture is much more difficult with the HIRS slit and so 
observation P2480101 consisted of three individual exposures. The A1032 line was used 
to check target alignment during the exposures as it is observed in all four of the FUSE 
telescope channels and has a high count rate. The LiFl channel was aligned throughout, 
with only a small loss of flux in the third exposure. LiF2 showed a loss of flux of around 10% 
compared to LiFl. A signiflcant number of counts were only detected from SiCl in the third 
exposure, but there are rapid variations in the count rate during this exposure and so these 
data are not useful. SiC2 showed a relatively uniform A1032 count rate during exposures 2 
and 3, but very few counts in exposure 1. 

The loss of flux in the SiC channels was anticipated and so a second observation was 
obtained with AG Dra in the LWRS aperture, but with the LiF focal plane assemblies moved 
such that the target would not lie in the aperture. The effective areas of the SiC channels 
are lower than those of the LiFs, while the degree of astigmatism is greater and so the count 
rate per unit area on the detector is much lower. This results in the count rate in the A 1032 
line lying at a safe level. No loss of flux is evident during the single SiC-only exposure. A 
weak signal from the O VI lines was detected in the LiF channels due to scattered light from 
AG Dra. 



5. The Spectrum and Line List 

Fig. 1 displays the spectrum of AG Dra over the complete FUSE wavelength range 
assembled from the LiF spectra of the 2000 March observation and the SiC spectra of the 
2001 April observation. The key features are (i) a uniform continuum from the white dwarf; 
(ii) narrow, nebular emission lines from highly-ionized species; and (iii) interstellar absorption 
fines. 

The present work will largely focus on the emission lines in the spectrum. We briefly 
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FLUX [ x10-'^ erg cm-^ s"' A"' ] 




Fig. 1.— FUSE spectra of AG Dra from SiClB (top), LiFlA (middle) and LiF2B (bottom). 
Spectra from SiClA (1080-1090 A) and SiC2B (1085-1095 A) have been used to complete 
the wavelength coverage. The shaded profiles in the middle plot show the Ovi and Hell 
hnes reduced by a factor of 50. 
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discuss some properties of the interstellar lines here. Comparisons of the absorption line 
wavelengths with those of the emission lines shows that they are redshifted relative to the 
radial velocity of the system by ~ 110 km s~^ and so are from the interstellar medium. 
Absorption by the hydrogen Lyman series, Ol, Nl, Nil, Fell, Sill and Arl is seen, as well 
as numerous lines of H2. Only one component to the absorption lines is seen; in particular 
no H2 absorption is seen at the radial velocity of the star. In the outer layers of the star 
one may expect the AG Dra nebula to cool sufficiently to form a significant layer of H2. 
We estimate an upper limit of 10^^ cm~^ to the H2 column density of this layer from our 
non-detection. 

The interstellar H2 absorption spectrum shows no lines from levels with J' > 5, indicat- 
ing relatively low excitation temperatures (< 1000 K). The H2 fines serve as useful fiducial 
marks for fixing the wavelength scale of the FUSE spectra as many of the fines are not 
saturated (unfike the atomic species). We adopt a velocity of —23 km for the interstellar 
lines based on lUE spectra (Viotti et al. 1983). 

In the following sections some references are made to the interstellar H2 lines, and 
we make use of the shorthand notation for referring to the H2 transitions introduced by 
Barnstedt et al. (2000). Thus, e.g., L13R3 refers to the u' — 13 to u" — 0, R(3) transition 
of the H2 Lyman band; while W0R3 refers to the u' — to u" — 0, R(3) transition of the 
Werner band. 



5.1. Line list 

The only AG Dra emission lines previously reported in the FUSE wavelength band have 
been the Ovi AA1032, 1038 fines, and Hell A1085 (Schmid et al. 1999). The complete list 
of lines found in the FUSE spectrum is given in Table 2. 

The measured wavelengths given for each line are those measured through Gaussian 
fits to the spectral lines, corrected to be placed in a reference frame where the interstellar 
absorption lines have a velocity of —23 km s~^ - see Sect. 12. The spectrum with the largest 
effective area at the particular wavelength considered was used to derive the wavelengths. 
The LiF measurements (990-1187 A) are all from the 2000 March observation, while the SiC 
measurements (905-1095 A) are from the 2001 April observation. Two fluxes are given for 
each line, corresponding to the two measurement dates. The first SiC measurements are 
more uncertain due to short exposure times. 

Line identifications have been based on solar spectra, in particular, the catalogs of 
Feldman et al. (1997) and Curdt et al. (2001). There are several lines in the AG Dra 
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spectrum for which no obvious candidates from the solar spectra could be found. Some of 
these lines are identified here and are marked with a (*) - further details can be found in 
Sect. 11. 



o 
O 




971.0 971.5 972.0 972.5 . 973.0 973.5 974.0 

Wavelength [ A ] 



Fig. 2. — SiC2A Spectrum from 2001 April showing the Nevil A973 emission hne. The 
spectrum has been binned by a factor 4. The shaded region denotes Ly7 airglow emission. 
The dashed and solid lines demonstrate that the emission line escapes absorption by the Lyy 
interstellar line (see text for details). Other interstellar absorption lines due to Ol and H2 
are indicated. 



6. Detection of Ne VII 

A weak hne is detected at 972.917 A which if corrected for the AG Dra systemic velocity 
using the measurement for the S VI A944 line (Sect. 12) gives a rest wavelength of 973.353 A. 
This matches the 2s 2p ^Pi - 2p^ transition of Nevil which has been detected in the 
laboratory by Edlen (1983) at 973.35 A and in solar spectra by Feldman et al. (1997). The 
transition is analogous to the Ov A1371 line which is formed, at the low electron tempera- 
tures found in the nebula, by recombination from O VI (Sect. 10). The identification of the 
Ne VII line thus indirectly implies the existence of Ne Vlll in the AG Dra nebula. 

The large radial velocity of AG Dra allows the Ne VII line to escape interstellar absorp- 
tion by the H2 LllPl transition at 973.348 A, yet it pushes the line closer to the interstellar 
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Table 1. Observing parameters for AG Dra observation 



Date 


Phase^ 


Datasct ID 


Exposure 


Aperture 


Time 


Exposure time (s) 


2000 March 16 


0.540 


S31202 


1 


LWRS 


15:57 


2388 


2001 April 25 


0.277 


P24801 


1 


HIRS 


05:28 


901 








2 


HIRS 


05:50 


788 








3 


HIRS 


06:09 


451 






P24802 


1 


LWRS 


06:57 


2060 



Using ephemeris of Fekel et al. (2000). 
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Table 2. AG Dra emission lines in the FUSE waveband. 



A ^ 

meas 


Flux'' (10- 


erg cm s ^ ) 




Ion 


Transition 


Wavelength 


V / 


Mar. 2000 


Apr. 2001 








Reference 


920.082 


- 


7 


920.561 


Hen 


Balmer (n=20) 


1 


922.336 


- 


7 


922.748 


Hen 


Balmer (n=18) 


1 


923.674 


- 


10 


924.147 


He II 


Balmer (n=17) 


1 


927.411 


18 


18 


927.851 


Hell 


Balmer (n=15) 


1 


929.886 




20 


930.342 


Hen 


Balmer (n=14) 


1 


932.989 


64 


48 


933.378 
933.449 


Svi 
Hen 


35 2 5i/2 - 3p 2p3/2 
Balmer (n=13) 


1 


936.929 


- 


11 


937.394 


Hen 


Balmer (n=12) 


1 


942.068 


35 


33 


942.513 


Hen 


Balmer (n=ll) 


1 


944.100 


19 


15 


944.523 


Svi 


3s 2Si/2 - 3p 2p3/2 




948.843 


- 


15 


949.329 


Hen 


Balmer (n=10) 


1 


958.241 


68 


56 


958.698 


Hen 


Balmer {n=9) 


1 


972.917 


- 


4 


973.35 


Ne vii 


2s 2p iPi - 2p2 i_D2 


3 


976.603 


9 


8 


977.020 


Cm 


2s^ - 2s 2p 


4 


989.388 


- 


9 


989.787 


Nin 


2s2 2p 2Pi/2 - 2s 2p2 ^0^12 


2 


991.078 


- 


7 


991.564 


Nin 


2s2 2p2p3/2 -2s 2p2 2£,g^2 


2 


991.795 


- 


24 


992.363 


Hen 


Balmer {n=7) 


1 


992.298 


18 


14 


992.731 


Nevi 


2s2 2p 2Pi/2 - 2s 2p2 4p3^2 


5 


996.671 


29 


26 


997.169 


Ne VI 


2s2 2p 2Pi/2 - 2s 2p2 4p^^2 


5 


998.812 


140 


135 


999.291 


Ne VI 


2s2 2p 2Pj/2 - 2s 2p2 4p^^2 


5 


1005.323= 


60 


52 


1005.789 


NeVi 


2s2 2p 2p3/2 - 2s 2p2 4p3^2 


5 


1009.899 


12 


12 


1010.323 


Nevi 


2s2 2p 2p3/2 _ 2s 2p2 4pj^2 


5 


1031.549 


5290 


3930 


1031.926 


Ovi 


2s 2Si/2 - 2p 2p,/2 


6 


1037.262 


1750 


1340 


1037.617 


Ovi 


2s 25i/2 - 2p 2p^/2 


6 


1062.213 


4 




1062.664 


Siv 


3s2 3p 2pj^2 - 3s 3p2 '^D.ij2 


2 


1065.707 


5 


6 


1066.195 


Fein 


3S> -^Gr, - 3(P {■^G)4p •'■F4, (*) 


7 


1072.527 


13 


6 


1072.973 


Siv 


3s2 3p 2P3/2 - 3s 3p2 2£)g^2 


2 


1073.068 


3 


5 


1073.518 


Siv 


3s2 3p 2P3/2 - 3s 3p2 ■^D3/2 


2 


1084.004 


240 


229 


1084.940 


Hen 


Balmer (n=5) 


1 


1124.264 


6 


8 


? 








1136.003 


27 


18 


1136.532 


Nev 


2s2 2p2 3p^ _ 2s 2p3 


8 


1136.494 


4 


3 


? 








1136.775 


4 


5 


? 








1138.971 


4 


2 


1139.494 


Fen 


(3p)4s b4P7/2 - 5p 4^5/2 (*) 


7 


1140.653 


8 


6 


1141.172 


Fen 


4s b4p5/2 - 5p 4Dg/2 (*) 


7 


1142.429 


9 


7 


1142.956 


Fein 


3^*^ 3£)3 _ (4(5)4p (^-j 


7 


1145.080 


62 


45 


1145.615 


Nev 


2s2 2p2 3p2 - 2s 2p3 


8 


1152.102 


8 


7 


? 









''Wavelengths are accurate to around ±0.02 A. 

''Absolute fluxes accurate to 15 % for strong lines, up to 40 % for weaJjest lines. 
"Wavelength was derived by estimating the center of the base of the line. 



Note. — Wavelength references; (1) Appendix A.3; (2) NIST atomic database; (3) Edlen (1983); (4) Boyce & Rieke (1935); 
(5) Appendix A.2; (6) Kaufman & Martin (1989); (7) Kurucz; (8) Appendix A.l. 
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H I Ly7 absorption. Fig. 2 shows a comparison of the SiC2A spectrum with a model spec- 
trum, demonstrating that the NeVll hne is not significantly absorbed by Ljj. The model 
shows a uniform continuum, and the NeVll A973.35 line, shifted by —133 km s~^. Two 
components were required to model the Ly7 absorption. The first is at a velocity of —23 km 
s~^ with a column density of 5 x 10^° cm~^ (a b value of 10 km s~^ was chosen). This was 
estimated by fitting the damping wings of the Hi Ly/3 absorption line. The second compo- 
nent is at a velocity of —108 km s~^ and has a column density of 5 x 10^^ cm~^. This is 
required to fit the absorption between 972.1 and 972.3 A and give the sharp edge at 972.1 A. 

No atomic data exist in the literature for interpreting the strength of the NcVII line. 
We note that the flux is around a factor 10 weaker than the Ov A1371 line (Mikolajewska 
et al. 1995). The solar photospheric 0/Ne ratio is 5.6 (Grevesse & Sauval 1998) and so the 
Nevil flux is consistent, within an order of magnitude, if the oxygen nebula conditions and 
atomic data are extended to neon. 

The imphed existence of Nevill in the AG Dra nebula is of great importance as the 
ionization potential of Ne+^ is 207.3 eV making it the highest ionization species ever found 
in the spectrum of AG Dra (e.g., the ionization potential of 0+^ is only 113.9 eV). It thus 
probes a much deeper part of the nebula. 

7. The O VI lines 

The two Ovi hnes are the strongest hnes in the FUSE spectrum and both show two 
distinctive features: asymmetric proflles and enhanced continuum at the bases of the lines. 
The A1038 line is affected by Cll A1037.018 and H2 L5R1 A1037.146 interstellar absorption 
on the short wavelength side of the proflle, which serves to make the A1032/A1038 flux ratio 
larger than the optically thin ratio of 2. The long wavelength side of A1038 is unaffected by 
the ISM and closely matches the shape and position of the A1032 line in velocity space. 

7.1. Line profile 

Fig 3(a) shows the A1032 hne plotted on a logarithmic flux scale, which reveals a P 
Cygni-hke proflle for the line. On this plot we have overplotted a Gaussian proflle derived 
from fltting only the continuum and the red side of the A1032 proflle. Also shown is the 
expected position of the L6P3 H2 absorption line. Clearly there is a possibility that it is the 
H2 line giving rise to a false P Cygni profile. Fig 3(b) shows the actual spectrum divided by 
the computed Gaussian spectrum, revealing the absorption to be ~ 200 km s~^ wide. This 
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ls much larger than the width of other H2 hnes in the spectrum - e.g., the nearby L6R3 
A1028.986 hne has a FWHM of 30 km s~^ - and so we attribute the absorption to the Ovi 
profile. The fact that the blue wing of the A1032 line is being absorbed (in addition to the 
continuum) shows that the O VI emitting region is an expanding, optically thick medium. 

A P Cygni profile in the N v A1238 hne seen in spectra from lUE was first reported by 
Viotti et al. (1983), and discussed later by Viotti et al. (1984), Lutz ct al. (1987) and Kafatos 
et al. (1993). In all these cases the P Cygni profile was only found during the outburst phase 
of AG Dra. This may simply be due to the low signal-to-noise in the continuum during 
quiescent periods, however, and Viotti et al. (1984) note that the A1238 line ". . .is always 
characterized by a sharper blue wing. . ." , an asymmetry consistent with the profile found 
for Ovi A1032 here. The significance of the 06 profile in relation to the structure of the 
nebula is discussed in Sect. 14. 




-300 -200 -100 100 200 300 -300 -200 -100 100 200 300 

Velocity [km/s] Velocity [km/s] 



Fig. 3.— (a) The LiFlA Ovi A1032 from the 2000 March observation, plotted on a loga- 
rithmic scale, reveahng the P Cygni absorption on the short wavelength side of the emission 
line. The expected location of the H2 L6P3 interstellar absorption is indicated. A Gaussian 
profile fitted to the continuum level and the long wavelength half of the profile is overplottcd. 
(b) The Gaussian profile shown in (a) is divided through by the observed profile, revealing 
the full extent of the absorption. 
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Fig. 4. — The LiFlA spectrum from 2000 March showing the broad wings of the O VI A1032 
and A1038 hnes. Overplotted is a model for the electron scattering wings assuming an 
electron temperature of 30,000 K and an optical depth of 0.04. 



7.2. Electron scattering wings of the O VI doublet 

The broad wings to the O VI profiles were measured previously by the ORFEUS echelle 
spectrograph (Schmid et al. 1999). They can be compared to the wings measured for the 
Hen A1640 line by lUE (Viotti et al. 1983). Possible instrumental causes of the wings (e.g., 
grating scattering or 'bleeding' of the detector electronics) can be ruled out by the presence 
of sharp interstellar absorption lines in the wings and the fact that the absorption in the C ll 
interstellar lines goes to zero. Schmid et al. (1999) suggested that the wings are caused by 
Thomson scattering of Ovi photons from nebular electrons, and equated the width of the 
wings to an electron temperature of 30 000 K. Using the expression of Bernat & Lambert 
(1978) for electron scattering in a layer outside of the line formation region, together with 
the Thomson scattering redistribution function of Mihalas (1970) we find an optical depth of 
0.04 and an electron temperature of 30 000 K fits the O VI AA1032, 1038 line wings reasonably 
well (Fig. 4). These values of the optical depth and electron temperature can be used to 
estimate a size and density of the O VI region. Since r aT^Ngr, where r is the optical depth, 
(Tt the Thomson cross-section, A^e the electron density and r the radius of the (assumed) 
spherical O VI emitting region, and the fiux in the Ovi A1032 line is related to and r as 
(assuming the line is optically thin) 
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And'F = 0.83'-^Ab{O)fA,,N!^ (1) 

where d is the distance to AG Dra, F the hne flux, the factor 0.83 represents the ratio of 
protons to electrons in a fully ionized plasma, h the Planck constant, c the speed of light, A 
the wavelength, Ab{X) the abundance of element X relative to hydrogen, Uj the population 
of the emitting level of the A 1032 transition relative to the ion population, and Aji is the 
radiative decay rate for the transition. The 0/H abundance is obtained from Smith et al. 
(1996), who measured the photospheric abundance in the giant, while the value of rij was 
obtained from the model of Ovi in the CHIANTI database (Young et al. 2003; Dere et al. 
1997), assuming an electron temperature of 30 000 K, and an electron density of lO^'^ cm~^ 
(note that nj/N^ is independent of A^e for the A1032 line, and so this choice does not affect 
the derived result). The Aji value was also obtained from CHIANTI. 

The two relations then imply 

N^r ^ 7.52 X 10^^ cm'^ (2) 
N^r^ ^ 4.15 X 10^^ cm-^ (3) 

Thus one has r 17.3 x 10^^ cm (= 11 Rq) and A^e ~ 1.0 x 10^^ cm~^. Given the uncertainties 
in the various parameters, we estimate that these quantities are accurate to no better than 
a factor of 5. The density is an order of magnitude higher than that derived from the O IV 
density diagnostic (App. B), which one may expect if the O VI is formed in a deeper layer of 
the nebula. 

This same calculation was performed by Viotti et al. (1983) when interpreting the 
broad wings of the Hell A1640 line observed by lUE during an active phase of the system. 
They ruled out the possibility of electron scattering on account of the large electron density 
(5.0 X 10^^ cm~^) they derived which was inconsistent with the density of ~ 10^°cm^^ derived 
from the lUE O IV density diagnostic (sec also Appendix B). We note that using the 2500 pc 
distance to AG Dra instead of their 700 pc distance increases their He ill region radius to 
7.7 X 10^^ cm and decreases their density to 3.9 x 10^^ cm~^, in reasonable agreement with 
the Ovi numbers. The optical depth of the electron layer is much greater from the Hell 
line than from the Ovi lines, although one may expect the opposite case as Ovi will be 
formed much deeper inside the nebula. However, Mikolajewska et al. (1995) found the AG 
Dra nebula to be larger during outburst, when the Viotti et al. (1983) measurement was 
made, which would thus lead to a larger optical depth. 

Another possibility for the broad wings is that they represent the thermal emission 
of a fast wind with speeds 1000-2000 km s~^. Such a scenario is presented in the AG 
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Peg observation of Nussbaumer et al. (1995), however in this case the broad wings are 
accompanied by a P Cygni profile in the Nv A1238 hne that extends out to velocities 
?a 1150 km s^^. For AG Dra, the P Cygni absorption reaches only ^ 200 km s~^, and so it 
is unlikely that the broad wings are due to a fast wind. In addition AG Peg has a history of 
high velocity Wolf-Rayet features (Kenyon et al. 1993) while AG Dra has no such history. 

8. Ne VI lines 

The five Ne VI 2s^2p ^Pj - 2s2p^ ^Pj' intercombination lines are prominent in the AG Dra 
spectrum between 992 and 1010 A. The A992 line is the weakest of the multiplet and found 
in the wing of Hell A992.363. Due to the large radial velocity of AG Dra the Ne VI A1005.696 
line is almost exactly coincident with the L8P2 H2 1005.397 A interstellar absorption line, 
leading to a double-peaked profile for the emission line. The Ne VI A1010.247 line is truncated 
on the short wavelength side of the profile by the H2 WOQl A1010.132 absorption. The A997 
and A999 lines are free from interstellar absorption and blending, and their ratio potentially 
allows a constraint on the nebular electron density to be made (Espey et al. 1996). 

The model we use for the Ne VI ion employs the electron excitation and radiative decay 
rates from v.4 of the CHIANTI atomic database (Young et al. 2003; Dere et al. 1997). 
Proton excitation rates are also included, but are negligible for the low temperatures found 
in nebular conditions. An important process to account for is photon excitation by the white 
dwarf radiation field. A star of effective temperature 1-2 x 10^ K emits predominantly in 
the extreme ultraviolet, between wavelengths 100-600 A. The strong 2s^2p-2s2p^ resonance 
transitions of Ne VI have energies between 400-600 A and so will be strongly pumped by such 
a radiation field at distances close to the star. 

Photon excitation can be accounted for in v.4 of CHIANTI through specifying a black 
body radiation field of temperature T^,, and the distance from the star's centre, r, measured in 
stellar radii units (Young et al. 2003). For the AG Dra model, we consider two temperatures: 
the value of 100,000 K derived by Mikolajewska et al. (1995) from fits to the UV continuum 
measured by lUE; and the value of 170,000 K derived by Greiner et al. (1997) from X-ray 
data from ROSAT. An electron temperature of 25,000 K is assumed. Figs. 5a,b show the 
effects on the Nevi A997/A999 ratio for each temperature. The four curves in each plot 
correspond to different values of r. The shaded region on each plot denotes the measured 
ratio, with la error bars. It is clear that the observed A997/A999 ratio is incompatible with 
the NeVI ions being located at <500 white dwarf radii (/?wd) from the white dwarf, and 
thus photoexcitation has a weak effect on the level populations. For a white dwarf radius of 
0.1 Rq (Mikolajewska et al. 1995), this places the NeVi ions at least 50 Rq from the white 
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dwarf and a significant way towards tfie giant at around 300 Rq (Mikolajewska et al. 1995). 
The ionization potential of Ne"*"^ is 126.2 eV, which compares with 113.9 eV for O"*"^, thus 
implying that the O VI emission is from regions even further from the white dwarf. One can 
thus rule out the suggestion (e.g. Viotti et al. 1983; Mikolajewska et al. 1995) that Hell and 
N V emission is from very close to the white dwarf. 

Our Ne VI model can be affected if the actual white dwarf spectrum deviates significantly 
from a black body. In particular, the very transitions that are pumped in the NeVI ion are 
likely to appear as absorption features in a hot white dwarf atmosphere. A Kurucz spectrum 
calculated for an effective temperature of 100,000 K and log g = 8.0, with solar abundances, 
shows Ne VI absorption lines whose depth at line center is around 30-60% of the continuum 
level. Thus photoexcitation may be overestimated by a factor of 2-3 in the black body 
model. We thus suggest that the A997/A999 ratio constrains the Nevi ions to be formed at 
a distance of at least 300 white dwarf radii. 

In the case of a neghgible radiation field, the electron density is restricted to the range 
5.5 < log A^e < 10.5. Increasing the radiation dilution factor narrows the allowed densities 
towards the higher end of this range, with the case of r = 500 implying densities of ~ 
10^° cm~^. This is consistent with densities derived from the intercombination transitions 
of OlV, measured with lUE Viotti et al. (1983). Fluxes obtained with lUE are reassessed 
in Appendix B with the most recent atomic data from CHIANTI. A density of 10^° cm~^ is 
derived. 



9. The Hell recombination lines 

The Balmer series of Hell (transitions n— >2, n>3) provide many emission lines in the 
AG Dra spectrum, from the n=5 line at 1085 A to the n=20 line at 920 A. The even members 
of the sequence lie close in wavelength to the Hi Lyman series, at around 0.3 A shortward of 
the Hi lines. This, combined with the ~ 100 km s^^ blueshift of the stellar lines compared 
to the interstellar absorption, is sufficient to make all lines from n=6 and above observable. 
Coincident blends with other interstellar absorption species prevent the strongest n=6 and 
n=8 lines being observed, but the n=10, 12 and 14 lines are all visible. These three lines 
are factors of 2-3 weaker than expected from the odd members of Balmer sequence and are 
blueshifted relative to these lines by Ki5-10 km s~^ (Sect. 12). This behaviour is puzzling. 
The reduced fluxes are not due to an intrinsic property of the He ll ions as recombination 
does not favor odd over even members of the sequence. Any absorption (either stellar or 
interstellar) must occur on the long wavelength sides of the lines to give rise to the observed 
blueshifts. The natural candidate is absorption through the damping wings of the hydrogen 
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Fig. 5. — Predicted variation of the Ne VI A997.14/A999.26 emission line ratio in the presence 
of 100,000 K (top panel) and 170,000 K (bottom panel) blackbody radiation fields. Four 
cases are shown, corresponding to the Nevi ions being different distances (r, in white dwarf 
radii units) above the siuface of the white dwarf. The grey band denotes the measured 
A997/A999 ratio from the 2000 March FUSE spectra. Atomic data are from the CHIANTI 
database. 
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lyman lines. However, both the A937 and A949 hnes show continuum on their long wavelength 
sides which is at the same level as the continuum in the rest of the spectrum, ruling out 
absorption by the H I lines. 

The strongest of the Hell hnes at A 1085 is affected by interstellar absorption in the Nil 
AA1084.566, 1084.584 lines which result in an asymmetric appearance to the A1085 profile. 
We estimate that the absorption removes around 30% of the He ll flux. 

The He ll Balmer lines can be used to place constraints on the density and temperature 
of the nebula. Storey & Hummer (1995) provide emissivities for all of the hnes up to 
n=50, calculated assuming Cases A and B of Baker & Menzel (1938) over a wide range of 
temperatures and densities. Case A requires the Hell Lya line to be optically thin, which 
is not applicable for the AG Dra nebula on account of the high densities. We consider 
the lines A958, A942 and A927 (n=9, 11, and 15, respectively) which are unblended and 
not contaminated by interstellar absorption. The closeness of the three lines in wavelength 
minimizes any distortions in the line ratios due to dust extinction. 

Fig. 7 plots the theoretical variation of the (A927+A942)/A958 ratio with density for 
a wide range of temperatures. Also plotted is the observed ratio derived from combining 
the SiClB and SiC2A spectra from the 2001 April observation, and which has a value of 
0.95 ± 0.08. It is clear that low nebula temperatures (T < 10^ K) are ruled out unless the 
density is also low (A'e < 10^ cm~^), while for hotter nebula temperatures low densities 
are ruled out. Other plasma diagnostics (Sects. 10 and Appendix B) suggest a density of 
10^° cm~^ and a temperature of 2 x 10^ K and these values are in reasonable agreement 
with the He ll lines. 

The Hell lines can also be useful in determining the extinction towards AG Dra. We 
consider the A958, A1085 and A1640 hnes (n=9, n=5 and n=3, respectively), the latter of 
which was observed many times with lUE and was found to have a flux between 2 x 10~^^ and 
8 X lO"^"*^ ergcm~^s^^ during the quiescent phase of AG Dra (Gonzalez- Ricstra et al. 1999). 
Assuming no extinction and nebula conditions of T = 2 x 10^ K and A^e = 10^° cm^^. Storey 
& Hummer (1995) predict A1640/A958 and A1640/A1085 ratios of 24.2 and 5.49, respectively 
Including the full range of variability of the A 1640 line gives upper and lower limits of 143 and 
36 for A1640/A958, and 24 and 6.0 for A1640/A1085 (adjusting the A1085 flux for intersteUar 
absorption). These values show the observed A958 and A1085 lines are weaker than expected 
from the A1640 line, consistent with a non-zero Eb-v- Assuming a more precise A1640 flux 
of 3.5 X 10~^^ based on the orbital phase of the 2001 April observation, and the extinction 
curve of Fitzpatrick (1999) allows us to estimate an E{B — V) value of 0.10 for AG Dra. 
This value is significantly higher than that found by Mikolajewska ct al. (1995), but we note 
that the Fitzpatrick (1999) extinction curve is an extrapolation to FUSE wavelengths based 
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on observations made above f=:il200A and so may be in error. 

10. The Nebula Temperature 

Methods for determining the electron temperature in the nebulae of symbiotic stars 
from emission line ratios are discussed in Sect. 3 of Nussbaumer (1987), with reference to 
lUE spectra. Temperatures for the nebulae were typically found to be ^20,000 K, indicating 
that the radiation field from the hot component of the symbiotic dominates the ionization 
balance of the plasma, rather than electron collisions. The FUSE wavelength range opens 
up further diagnostic possibilities, and the use of O VI and C ill lines to constrain the plasma 
temperature are discussed below. 

10.1. O VI, O V temperature diagnostic 

In a highly-ionized photoionized plasma such as the AG Dra nebula, it is possible to 
use the Ovi A1032 and O v A1371 emission lines to determine the nebula temperature. For 
electron temperatures ~ 10^ K the 2s2p ^Pi - transition of beryllium-like ions (e.g., 

O V A1371) is predominantly formed through dielectronic recombination from the lithium-like 
ion (e.g., Ovi) into the level. The cmissivitics of such recombination lines vary slowly 
with temperature. In contrast, for emission lines excited directly through electron collisions, 
the emissivity falls sharply with temperature on account of the MaxweUian cxp(— AE/ZcT) 
term in the expression for the excitation rate coefficient. For temperatures > 10^ K, the 
A 1032 line is sufficiently strong that it remains electron-excited rather than recombination- 
excited, and so the Ovi A1032/OV A1371 ratio is an excellent temperature diagnostic, 
independent of the relative concentrations of the two ions. 

Recombination data suitable for predicting the strength of the O V A1371 line are from 
Nussbaumer & Storey (1984), while data for the A1032 fine are from v.4 of the CHIANTI 
database (Young et al. 2003). The A1032/A1371 ratio varies by several orders of magnitude 
over a small temperature range (Fig. 8), providing a tight constraint to the electron temper- 
ature. Measurements of the A1371 line at different orbital phases of AG Dra were presented 
by Mikolajewska et al. (1995). Using these values together with the Ovi A1032 flux from 
FUSE gives the ratio indicated in Fig. 8, yielding a temperature of ~ 25, 000 K, in good 
agreement with the temperature of 30,000 K derived from the electron scattering wings on 
O VI fine profiles. 
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10.2. cm 

Cm gives rise to two dominant spectral features in the FUSE waveband: the 2s - 
2p transition at 977.020 A and the set of six 2s2p ^Pj - 2p'^ ^Pji transitions between 1174 
and 1176 A. These hnes together with the A1909 line observed by lUE potentially allow both 
the temperature and density of the plasma to be estimated. 

The A977 profile from the 2001 April observation is displayed in Fig. 9 and shows 
the emission line to lie within a deep absorption trough in the white dwarf continuum. The 
absorption on the long wavelength side of the profile is due to C ill in the interstellar medium 
- the AG Dra emission fine is saved on account of the high radial velocity of the system. To 
the short wavelength side is the Ol A976.448 interstellar absorption fine. The centroid of 
the emission line is redshifted by 5-15 km s~^ relative to other species indicating that it 
is partially absorbed by the ISM components. 

The A1176 feature is not detected (Fig. 9) and we estimate an upper fimit to the flux 
of 4 X IQ-^'' ergcm-^s-^ The A1909 was detected by lUE, and Mikolajewska et al. (1995) 
give quiescent fluxes of 1-4x10"^^ ergcm~^s~^. 

Fig. 10 shows the theoretical variation of the A1176/A977 and A1176/A1909 ratios, 
derived from the CHIANTI database. In each case the ratios are plotted at two temperatures: 
80,000 K and 25,000 K, the former corresponding to the temperature of maximum abundance 
of Cm assuming an electron-ionized plasma, and the latter the temperature derived in 
Sect. 10.1 for the O VI region of the AG Dra nebula. 

The measured A977 flux is a lower fimit as the extent of interstellar absorption is not 
known. The resulting A1176/A977 ratio can thus not discriminate the temperature of the 
plasma. Using the lowest of the measured A1909 fluxes, the upper limit to the A1176/A1909 
ratio is overplotted in Fig. 10(b). Clearly the ratio is incompatible with a temperature 
of 80,000 K. A deeper exposure of the A1176 line together with a simultaneous and more 
accurate measurement of the A1909 line are required to improve the constraints on nebula 
temperature. 

11. Fluoresced and unidentified lines 

There are several emission lines in the AG Dra spectrum for which no obvious identifica- 
tions could be found based on solar spectra. These lines are listed in Table 3, together with 
their wavelengths corrected by a velocity shift of +136 km s~^, which is the velocity of the 
NeV A1136 line. The corrected wavelengths should be within ±10 km s^^ of the laboratory 
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wavelengths of the transitions. 

As many unidentified emission lines exist in the UV spectrum of the symbiotic nova RR 
Telescopii that were later identified as Fe ll transitions fluoresced by the strong nebular hues 
(Johansson 1988), we searched for Fell lines in the FUSE bandpass that could be similarly 
excited. We initially searched for the lines identified by Harper et al. (2001) who found many 
Fell lines fluoresced by H I Lya in the FUSE spectrum of the hybrid supergiant a TrA. None 
of the imidentified FUSE lines matched the Harper et al. list. In particular, the 1133.70 A 
and 1139.02 A lines, pumped via a line only 0.182 A from Lya line center, are not present 
implying either a narrow Lya line of FWHM ^ 0.36 A, or that the Lyo; line is formed far 
from the surface of the giant. 

In order to find other fluorescence lines in the FUSE waveband a simple model was de- 
veloped to predict spectra arising through radiative pumping of Fe II and Fe III by a Gaussian 
shaped emission line at a specific wavelength. The model calculates the quantity 

(4) 

where i, j and k are indices for the Fe II levels such that the transition i k is the transition 
pumped by the emission line and k ^ j is the subsequent radiative decay; 0(A) is the 
intensity of the fluorescing emission line proflle at wavelength A; Ni is the population of level 
i relative to the ion population; A^j the radiative decay rate for transition k j; and uJi is 
the statistical weight of level i. The model allows the centroid and width of the fluorescing 
emission line to be varied. The are distributed according to the Boltzmann distribution, 
and the energy level and radiative data for Fell and Felll are from the compilations of 
R.L. Kurucz, available at http://cfaku5.harvard.edu/. 

We considered pumping of both Fell and Felll by the O VI AA1032, 1038 doublet, and 
Table 4 shows the model predictions. A A1038/A1032 ratio of 1:2 was assumed, and the full 
width at half maxima of the O VI lines was taken as 0.30 A. For both Fell and Felll there is 
only one significant transition that is pumped by the O VI lines; the resulting fiuoresced lines 
thus have relative fluxes determined purely by the radiative decay rates for the transitions. 

The Fell hues fluoresced by A1032 are all pumped through the same Fell transition. 
The A1776, A1881 and A1884 hues have been identifled by Johansson (1988) in lUE spectra 
of RR Tel and have also been measured in lUE spectra of AG Dra by Viotti et al. (1983) 
although the authors did not identify the transitions. The strongest of the predicted lines 
below 1200 A is at 1141.172 A which provides an excellent match with one of the AG Dra lines 
(Table 3). Another good wavelength match is found for the A1139.494 transition, however. 
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Fig. 6. — Plots showing the Hell A937, A949 and A1085 emission hnes (Balmer n=12, 10 
and 5, respectively). Figures (a) and (b) show continuum emission on the long wavelength 
side of the emission lines, demonstrating that the reduced fluxes in the lines are not due to 
interstellar hydrogen absorption. Figure (c) shows the expected positions of N ll transitions, 
both in the cases of them being due to interstellar absorption or being in emission from the 
star. 



Table 3. Wavelengths for 
unidentified lines. 



Measured 


Corrected 


wavelengths 

(A) 


wavelengths 

(A) 


1065.752 


1066.141 


1136.494 


1137.010 


1136.775 


1137.291 


1138.971 


1139.488 


1140.653 


1141.171 


1142.429 


1142.948 


1152.102 


1152.625 
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Fig. 7. — The lines show the predicted variation of the Hell (A927+A942)/A958 ratio with 
density for a range of electron temperatures, derived from the calculations of Storey & 
Hummer (1995). The shaded region shows the ilcr errors on the observed ratio. 

the predicted A1139/A1141 ratio is much lower than the observed ratio. None of the other 
transitions can be identified in the AG Dra spectrum. 

Applying the same model to Fe III yields two lines whose wavelengths are close to mea- 
sured lines: A1066.195 and A1142.956. The latter is the strongest of the predicted lines in 
the FUSE waveband and is in excellent wavelength agreement with the line measured at 
1142.429 A. The predicted A1070.557 line would be expected at a wavelength of ps1070.16 A 
which is coincident with the H2 L3P3 A1070.142 interstellar absorption hne and so is not 
seen. Although the model A1066.195 line is expected to be much weaker than the A1070.557 
and A1142.956 hnes, there is a weak line found at a corrected wavelength of 1066.166 A which 
we identify with this transition. 

We conclude that two of the unidentified hnes can be identified as Fell A1141.172 and 
Felll Al 142.429, fluoresced by Ovi A1032. Two more lines can be tentatively identified as 
Fell A1139.494 and Felll A1066.195 based on wavelength matches, although the theoretically 
predicted fluxes are weaker than the observed values. As Fell and Felll are extremely 
complex ions, there is the possibility that the atomic data used in the theoretical model may 
cause the discrepancy in the fluxes. 
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Table 4. Predicted fluorescence lines. 



Pumping line Pumped transition Fluoresced transitions Relative flux 

(A) (A) 



Fluoresced Fell lines 



OviA1032 1032.041 



OviA1037 1037.808 



978.697 


6.8 


981.468 


5.5 


1114.355 


1.8 


1126.677 


2.8 


1139.494 


1.3 


1141.172 


10.0 (*) 


1142.370 


2.8 


1776.421 


2.5 


1776.660 


14.4 


1881.201 


16.7 


1884.116 


11.3 


1448.393 


1.3 


1663.782 


0.7 



Fluoresced Felll lines 



Ovi A1032 


1032.119 


1066.195 


1.5 






1070.557 


9.1 






1142.956 


10.0 ( 






2097.696 


30.4 


Ovi A1037 


1037.459 


1035.769 


1.4 






1071.747 


0.6 






1142.461 


0.7 






2103.809 


1.9 
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Fig. 8. — Theoretical variation of the O VI A1032/O V A1371 ratio. The dashed hnes indicate 
upper and lower hmits to the observed ratio based on the present FUSE observations of the 
A1032 line and lUE observations of the A1371 line. 
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Fig. 9. — SiC2A spectrum from 2001 April, showing the Cm A977 emission line. The 
absorption features on each side of the emission line are due to interstellar absorption by O I 
A976.448 and Cm A977.020. 
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Fig. 10.— Plots of the Cm A1176/A977 and A1176/A1909 ratios as a function of density 
at temperatures of 25 000 K (upper curve) and 80 000 K (lower curve). The horizontal lines 
show upper hmits derived from FUSE and lUE spectra. Calculated using the CHIANTI 
atomic database. 
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12. Emission line widths and shifts 

Fricdjung et al. (1983) reported systematic redshifts of resonance lines relative to in- 
tercombination lines from lUE spectra of seven symbiotic stars (including AG Dra), with 
magnitudes ranging from 10 to 40 km s~^. The values given for AG Dra were 12 to 15 km 
s~^. Viotti et al. (1984) gave a redshift of 9 km from their quiescent lUE spectra, but 
also noted that the N v AA1238, 1242 doublet exhibited a larger redshift of 30 km s^^ relative 
to the intercombination lines. The first direct measurement of the O VI lines by Schmid et 
al. (1999) yielded a redshift of 49 km s~^ relative to the system radial velocity for the A1032 
line. 

The cause of these redshifts has been attributed to the scattering of hne photons in an 
expanding medium (Priedjung et al. 1983), and this also explains the asymmetries noted in 
lUE spectra of the N v lines (Mikolajewska et al. 1995) and the O VI lines found here. 

Before discussing the FUSE results, the accuracy of the FUSE wavelength scale must 
first be considered. Due to the lack of an on board calibration lamp, the absolute wavelength 
cahbration of FUSE is good to < 70 km s~^ for observations with the LWRS aperture due 
to motions of the target within the aperture, and the thermally-induced grating rotation^. 
The relative wavelength calibration within a single spectrum is accurate to around 5-10 km 
s~^. In order to determine the absolute wavelength scale of the spectra it is thus necessary to 
measure features in the spectra for which the wavelengths are known through other methods. 
Typically one uses interstellar absorption lines for which velocities are known from lUE or 
HST spectra. For AG Dra we use the interstellar velocity measured by Viotti et al. (1983) 
of —23 km s~^ from lUE spectra using low ionization species. For each of the FUSE spectra 
several interstellar absorption lines are selected and their centroids are measured through 
line fitting. If the average velocity of these lines relative to their rest wavelengths is v, then 
we apply a correction of 23 — f km s^^ to that spectrum. The lines selected in each spectrum 
were based on the spectra of Barnstedt et al. (2000). 

The velocities of selected emission lines are presented in Table 5. For the asymmetric 
A 1032 line the centroid was estimated through minimizing the quantity 

^(A) = |A-A,|xF(A,) (5) 

where F{Xj) is the flux in the wavelength bin Xj. 

The Nev and Nevi lines are intercombination transitions, and are found at velocities 



^The FUSE Wavelength Calibration White Paper, http://fuse.pha.jhu.edu/analysis/calfuse_wpl.htnil. 
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Table 5. Velocity shifts of FUSE lines. 



Ion 


Line 

(A) 


Velocity shift 
(km s~^) 


Comment" 


Hen 


920.561 


-155.7 


L18P2 A920.242 (r) 




922.748 


-133.9 






923.796 


-153.4 






927.851 


-145.4 






930.342 


-146.9 






937.394 


-148.7 


L16P3 A936.859 (b) 




942.513 


-141.5 






949.329 


-153.5 






958.698 


-142.9 




b VI 


944.523 


— 134.3 




Cm 


977.020 


-127.9 


Oi A976.448 (b), Cm A977.020 (r) 


Nevi 


992.731 


-130.8 






997.169 


-140.4 






999.291 


-143.7 






1005.789 


-138.9 


L8P2 A1005.397 


Ovi 


1031.926 


-109.5 






1037.617 


-102.6 


Cii A1037.020 (b) 


Siv 


1062.664 


-127.2 






1072.973 


-124.6 




Nev 


1136.519 


-136.1 






1145.607 


-137.9 





^Interstellar absorption hnes that affect the velocity shift determi- 
nation are identified here. The letter in bracket denotes whether the 
absorption is to the blue (b) or red (r) side of the stellar emission line. 
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of —125 — 135 km s~^. Aside from the Ovi lines that are affected by P Cygni profiles, the 
other resonance lines are from S VI and S IV, which are consistent with the intercombination 
lines. Thus there is no indication from the FUSE spectra that the resonance lines are 
redshifted relative to the intercombination lines. The consistency of the S VI velocity shift 
with the intercombination lines reveals that this high ionization state is not showing a P 
Cygni absorption, restricting the size of the region participating in the wind. 

The Hell Balmer lines are all significantly blueshifted relative to the other emission 
lines, which may indicate that the model assumed in deriving the rest wavelengths is in error 
(Appendix A). 

The centroid of the O VI A1032 line is redshifted relative to the systemic velocity by 
39 km s~^ in reasonable agreement with the value of 49 km s~^ found by Schmid et al. 
(1999). 

Accurate inoasurcmcnts of emission line widths require good signal-to-noise relative to 
the continuum level, and so we give values only for the strongest lines in the spectrum in 
Table 6. The Ovi width represents that of the observed profile. One can estimate the true 
width of the line by assuming that it would have a shift comparable to the other emission 
lines in the spectrum of ~ —130 km s~^, leading to a width of ~ 105 km s~^. One may 
speculate that the broader width found for Ne VI compared to Ne V may be due to the Ne VI 
ions participating in the wind seen in the O VI profile. 



13. Time variability 

Including the present observations, the 900-1200 A region of the AG Dra spectrum has 
now been observed four times, and the fluxes of the prominent emission lines are given in 
Table 7. The ORFEUS-I/BEFS fluxes have been derived from the archived spectra available 

Table 6. Widths of selected FUSE emission hnes. 



Ion Wavelength (A) Width (A) Width (km s"^) 

Hen 958 0.177 55.4 

Nevi 999 0.257 77.1 

Ovi 1032 0.218 63.3 

Nev 1136 0.192 50.7 
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via the Multimission Archive at STScI (MAST)2, while the ORFEUS-II/TUES fluxes are 
taken from Schmid et al. (1999). 

Monitoring of AG Dra with lUE showed that the emission hnes are modulated by a 
factor two during the orbit, with the strongest fluxes at phase 0.5 ( Gonzalez- Riestr a et al. 
1999), when the white dwarf is in front of the giant. Both ORFEUS measurements of the 

Ovi and He II lines lie above the FUSE values, even though the 2000 March data were 
obtained close to phase 0.5. This suggests that the two ORFEUS measurements were not 
taken during periods of quiescence, and this is confirmed by inspection of Fig. 7 of Galis et 
al. (1999) where the U band flux is enhanced above quiescent levels at Julian dates 2449 249 
and 2450 409, corresponding to the ORFEUS observations. 

Comparing the line fluxes from the two FUSE observations shown in Table 2 shows that 
the Ovi lines show the largest decrease in the 2001 April data-set, being reduced by 25%. 
This is consistent with the factor of 2 modulation during the entire AG Dra orbit found by 
Gonzalez- Riestr a et al. (1999). The Ne VI hnes which have a comparable excitation potential 
to the O VI lines show a much smaller decrease in flux of around 10%. This may indicate that 
the changes in the O VI flux may signal variations of the optical depth with orbital phase. 

There was no variability beyond 2a levels observed during 2 ksec of the FUSE observa- 
tions of the Ovi hnes, despite there being 35 counts s~^ in the LiFlA channel. 

14. Discussion 

The previous sections have presented FUSE spectra obtained in 2000 March and 2001 April 
during a quiescent phase of AG Dra's ^15 year outburst cycle. These represent the first high 
resolution, high sensitivity spectra obtained in the 905-1 187 A spectral region for AG Dra. 
The high ionization states found in the FUSE bandpass reveal significant new information 
about the system, unavailable from earlier, longer wavelength lUE spectra. In particular, 
the presence of the Nevil A973 recombination line implies the existence of Nevill in the 
nebula, the highest ionization species ever recorded for AG Dra. The O VI A1032 line has 
a P Cygni profile, demonstrating a wind in the high ionization part of the nebula during 
quiescence. Previously P Cygni profiles had only been seen in lUE spectra obtained during 
outbursts. The NeVI intercombination lines are sensitive to the radiation field of the white 
dwarf and their ratios imply that the Ne VI emitting region (and thus also the O VI emitting 
region) is at least 300 white dwarf radii (or 0.14 AU for R^^ ~ 0.10 Rq) from the white 



^http://archive. stsci.edu/mast. html 
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dwarf. Evidence for the O VI being formed far from the white dwarf is also found from the 
Fe II and Fe ill fluorescence lines: the O VI resonance lines are capable of fluorescing iron 
lines in the giant's atmosphere, whereas Hi Lya is not, suggesting Ovi is formed close to 
the giant's surface, while Lya is formed much further out. 

The classical view of symbiotic nebulae has material flowing out from the giant and 
being ionized in the vicinity of the white dwarf. High ionization species (such as He ll and 
Ovi) are found in the region around the white dwarf (e.g. Friedjung et al. 1983), while 
low ionization species are formed at increasingly further distances. An alternative theory 
is that the strong nebular lines in high density systems such as AG Dra arc formed in the 
illuminated part of the giant's atmosphere, helping to explain the orbital modulation of the 
UV line fluxes found from lUE (e.g., Altamore et al. 1981). Proga et al. (1996, 1998) have 
created non-LTE photoionization models for the ionization of a giant's atmosphere in the 
case where it is illuminated by a hot companion star and have demonstrated that the general 
features of symbiotic star spectra can be reproduced with high effective temperatures and 
luminosities of the companion. The large emission line fluxes of highly-ionized species can 
only be reproduced if the giant has a wind extending to 2-3 giant radii. 



0.6 




0.01 0.10 1.00 10.00 

Distance from white dwarf [ AU ] 

Fig. 11. — Predictions for the location of the O Vl, Nevi, and Nevill ions in relation to the 
white dwarf, assuming a simple model for the ionization structure of the nebula (App. C), 
and a white dwarf black body temperature of 170,000 K. 

Evidence for the Ovi lines being formed in the giant's atmosphere comes from the 
aforementioned Nevi result and Fe fluorescence lines. We can also consider the location of 
the Nevi, Vlll and O VI ions in the system through ionization balance arguments. The high 
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Fig. 12. — A schematic diagram illustrating the locations of the O VI and NeVIII emitting 
regions (vertical lines) suggested by the present work. 

ionization potential of these ions means that they are produced through irradiation from 
the X-ray section of the white dwarf's radiation field. Fits to the measured X-ray spectrum 
(Greiner et al. 1997) are consistent with a blackbody field of temperature 170,000 K. By 
constructing a simple model of the ionization/recombination processes in the plasma (Ap- 
pendix C) we can calculate the fractional population of Nevi and Ovi as a function of 
distance from the white dwarf (Fig. 11). The distance here corresponds to the inner face of 
a nebular shell, and the model assumes that the radiation field is unattcnuatcd between the 
white dwarf and inner face. It is clear that such a hot blackbody prevents Ovi and NcVI 
being formed out to distances of 0.5-1.0 AU. For the model of AG Dra presented in Fig. 9 of 
Mikolajewska et al. (1995) this places the O VI and Nevi regions close to the surface of the 
giant, and well away from the white dwarf. Nevill, however, must be formed much closer 
to the white dwarf. A cartoon illustrating the structure of the AG Dra nebula is shown in 
Fig. 12. A difficulty for this particular model lies in the P Cygni profile of the Ovi line, 
since the continuum in the far ultraviolet is due to the white dwarf and not the giant: for 
white dwarf continuum photons to be absorbed by Ovi at orbital phase 0.5 (white dwarf 
in front of the giant), the Ovi ions are required to lie in front of the white dwarf along the 
observer's line of sight. We note however, that the electron scattered wings around the O VI 
hues significantly enhance the continuum (Fig. 4), and that the P Cygni absorption extends 
down to around 3 x 10~^^ ergcm~^s~^, the level of the white dwarf continuum, and not 
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to zero flux. We thus suggest that the O VI ions are absorbing the electron-scattered O VI 
photons and not the steUar continuum. This can be understood if the electron scattering is 
occurring in the densest part of the Ovi region close to the giant, with the absorption of 
the scattering wings then occuring in the lower density, extended wind. The high density of 
10^^ cm~^ derived from the analysis of the wings in Sect. 7.2 supports this idea. 

A full understanding of the structure of the AG Dra nebula can only come through 
detailed modelling of the ionization and wind structure, constrained by high quality spectral 
data. 

We thank S.D. Friedman for scheduling the 2000 March observation. P.R. Young thanks 
W.M. Moos for the allocation of observing time for the P248 program. J. Aufdenberg and 
R.L. Kurucz are thanked for useful discussions. 

A. Rest wavelengths of FUSE emission lines 
A.l. NeV 

The two NeV lines seen in the FUSE spectrum arise from decays of the 2p level to 
the ground 2s ^-Pi,2 levels. The separation of the two ground levels is strongly constrained 
by the measurement of the ^Pi - transition wavelength from SWS spectra obtained with 
ISO by Feuchtgruber et al. (1997). They find a wavelength of 14.3217±0.0002//m. 

The UV lines have only previously been measured by space-borne, solar UV spectrom- 
eters. Sandlin et al. (1977) give wavelengths of 1136.51±0.02 A and 1145.61±0.02 A. Edlen 
(1985) made use of these wavelengths to refine the energy levels of Nev, leading to revised 
wavelengths of 1136.51 A and 1145.59 A. More recently, the SUMER instrument on board 
SOHO has measured the lines. Feldman et al. (1997) give wavelengths of 1136.52 A and 
1145.62 A from off-limb spectra, while Curdt et al. (2001) give wavelengths of 1136.56 A and 
1145.66 A from disk spectra. The off-limb spectra are hampered by blending with coronal 
lines, though, while disk spectra of solar transition region lines typically show redshifts of 
5-10 km (e.g., Peter & Judge 1999). A detailed reassessment of SUMER disk spec- 
tra (Curdt 2001, private communication) yields revised wavelengths of 1136.551±0.020 A 
and 1145.632±0.010A. We correct these assuming a redshift of 5 km s ^ and use the error 
estimates from the UV and IR wavelength measurements to minimize the energy level sepa- 
rations of the three levels. These yield new wavelengths for the UV lines of 1136.532 A and 
1145.615 A which are used in the AG Dra analysis. 
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A.2. NeVI 

As with the NeV hnes, the NeVI intercombination hnes have only been measured in 
solar spectra. The most recent determination of the lines' wavelengths is that of Dwivedi 
et al. (1999) from off-limb SUMER spectra. However, these wavelengths show a serious 
discrepancy with the measurement of the Ne VI ground transition wavelength of Feuchtgruber 
et al. (1997) from the SWS instrument on ISO. In particular the Dwivedi et al. wavelengths 
for the - ^-Pi/2,3/2 transitions imply, when combined with the Feuchtgruber et al. (1997) 
wavelength, that the ^^3/2 - ^-Pi/2,3/2 transitions occur at 1005.777 A and 1010.303 A. Dwivedi 
et al. actually measure the latter lines at wavelengths 1005.696 and 1010.247 A - differences 
of 0.081 A and 0.056 A that are much larger than the 0.015 A accuracy quoted by Dwivedi 
et al. (1999). 

Analysis of more recent SUMER spectra (W. Curdt, private communication, 2001) with 
updated calibration files has resulted in revised wavelengths for the NeVI lines giving values 
of 992.731, 997.169, 999.291, 1005.789 and 1010.323 A. These values have been used in the 
AG Dra analysis. 

A.3. Hell 

Each of the Hell Balmer series lines consists of seven individual transitions that, in 
terms of wavelength, can be split into two groups. E.g., for the n=5 transitions the groups 
are at average wavelengths of 1084.912 A and 1084.977 A (based on data from the NIST 
database) - i.e., a velocity separation of 18.0 km s~^. In order to compare the velocity shifts 
of the He 11 lines with those from other ions, it is necessary to model how the individual line 
components contribute to the total line flux. 

Clegg et al. (1999) provides such a model for the Hell A1640 (n=3) and A1215 (n=4) 
lines, yielding estimates of the line centroids for a range of temperatures and densities, and 
for Cases A and B of Baker & Menzel (1938). We extend this model for 5< n< 20 with the 
method outhned in Sect. 3 of Clegg et al. (1999), atomic data from Storey & Hummer (1995), 
and ^4- values calculated with the FORTRAN routine described by Storey & Hummer (1991). 
Energy level data were obtained from NIST. For n>12, energy values were not available for 
all of the ^P, ^D and ^S levels. However, the wavelength splittings of the individual line 
components are largely determined by the splitting of the n=2 levels for these high n levels. 

To determine the expected wavelengths of the He ll Balmer hnes for AG Dra, an electron 
density of 10^*^ cm~^ and a temperature of 2 x 10^ K were assumed together with Case B. 
The emissivities and wavelengths of the seven individual line components were calculated. 
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and a combined synthetic profile was computed by assuming lines broadened by 50 km s~^. 
The resulting profile was then fit with a Gaussian, the centroid of this fit is the wavelength 
that is listed in Table 2. 



B. Revised O IV density from lUE data 

OlV belongs to the same isoelectronic sequence as Nevi, and the intercombination 
transitions that occur for Ne VI in the FUSE waveband are found for O IV between 1397 and 
1408 A. These lines form ratios that are sensitive to the electron density between 10^ and 
10^^ cm^"^, and Fig. 13 shows the ratios, relative to the strongest A1401 line, calculated from 
the CHIANTI database at an electron temperature of 25 000 K. Measurements of the O IV 
lines were made by lUE and the high resolution mode of the satellite allowed the individual 
transitions to be measured. Two data-sets with good measurements of the lines are from 
1981 August 3 and 1983 June 7. The former was obtained several months after the system 
went into a major outburst, while the latter was obtained during a quiescent phase. The 
line fluxes were very similar in each case, and the line ratios, relative to A1401, are shown 
in Table 8 and overplotted on Fig. 13. The ratios clearly show that the density lies between 
10^ and 10^° cm~^. The stronger lines have a higher signal-to-noise and so we favor densities 
closer to 10^° cm~^ and this is the value assumed for the density of the nebula in the present 
paper. 



C. Ionization model 

In Sect. 14 results are presented from a model of the ionization balance of oxygen and 
neon ions in a photoionized plasma. The details of this model are given here. 

The distribution of element ionization states in a plasma photoionized by a distant 
radiation source are determined by the distance to the source and the local density and 
electron temperature of the plasma. The number of ions leaving an ionization state i through 
photoionization is given by 

ai^ATrNiW{r) J B{n,X)a{X)dX (CI) 

where a is the photoionization cross-section, B is the specific intensity of the radiation field, 
and W{r) is the dilution factor given by 
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Table 7. Emission line fluxes from ORFEUS. FUSE 





Flux 1 


[x 10 erg cm 


-2s-l) 




ORFEUS-I 


ORFEUS-II 


FUSE 


Phase: 


0.723 


0.833 


0.036 0.277 


Line 








Ovi A1032 


8610 


6210 


5290 3930 


Hell A1085 


536 


248 


240 229 


Nevi A999 


448 




140 135 


Nev A1136 


34 




27 18 



Table 8. lUE O IV line ratios (energy units) 





SWP14641 


SWP20162 


Ratio 


1981 August 


1983 June 


A1399/A1401 


0.19 


0.25 


A1404/A1401 


0.34 


0.33 


A1407/A1401 


0.18 


0.19 
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Fig. 13.— Theoretical OlV line ratios, relative to the A1401.171 line, from the CHIANTI 
database, calculated for an electron temperature of 25 000 K. The symbols indicate the ratios 
from the 1981 August and 1983 June lUE observations, and the corresponding densities. 



W = 



1/2- 



(C2) 



and r = R^jR is the ratio of the radiation source's radius to the distance of the plasma from 
the source's center. 



The number of ions entering into the state i through recombinations from the i + 1 state 



IS 



N,i5, = N,^^Nla,J:r) + a^lT)) (C3) 

The set of linear equations 



N^{ai + A^eA-i) = ^i-iOLi-x + Ni+iN^Pi (C4) 
are solved to yield ion populations as a function of radius 
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Photoionization cross-sections for both the oxygen and neon ions are from the ana- 
lytic fits of Verner & Yakovlev (1995). The combined radiative and dielectronic recombi- 
nation rates for the oxygen ions are from Nahar (1999), while for the neon ions the ra- 
diative recombination rates arc from the fits provided by D.A. Verner and available from 
http://www.pa.uky.edu/~verner/fortran.html, while the dielectronic recombination rates 
are from Mazzotta et al. (1998). 
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